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0\ ■ Abstract 

, ^ , We study fluctuations of electric current in a quantum resistor and derive 

a general quantum-mechanical formula for the distribution of transmitted 
charge. For that we introduce a scheme of current measurement that in- 
volves a spin 1/2 coupled to the current so that it precesses at the rate pro- 
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T-|- ■ portional to the current. Our approach allows the study of charge transfer 

^P , without breaking the circuit. We analyze a single channel conductor and de- 

rive electron counting statistics for arbitrary relation between temperature 
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^D ' and voltage. For a perfectly transmitting channel the counting distribution is 






gaussian, both for zero-point fluctuations and at flnite temperature. At con- 
stant voltage and low temperature the distribution is binomial, i.e., it arises 
53 . from Bernoulli statistics. 
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Introduction. Instantaneous measurement is described in quantum mechanics by 
wavepacket reduction, i.e., by projection on eigenstates of an observable. A different kind 
of measurement, extended in tlie time domain, is realized in detectors and other counting 
devices. It is known that in such cases a certain revision of the measurement description 
is necessary. A famous example from quantum optics is the theory of photon detectors [|1| . 
Due to Bose statistics, the photons entering a photo-counter are correlated in time, and this 
makes the theory of photon detection a problem of many-particle statistics. For a single 
normal mode of radiation field the probability Pm to count m photons over time t is given 
by 

P™ = ^(:(a+are-''*'^^'^:), (1) 

ml 

where a+ and a are Bose operators of the mode, r] is "efficiency" of the counter, (...) stands 
for the average over a quantum state. The normal ordering : . . . : is important — physically, 
it means that, after having been detected, each photon is destroyed, e.g., it is absorbed in 
the detector. Instead of the probabilities, it is more convenient to deal with the generating 
function 

X(A) = Y. Pme'^"" ■ (2) 

m 

For the single normal mode Eq.(|l]) leads to 

X(A) = ( : expr/t(e^^-l)a+a: ) . (3) 

Eqs.([l|,0,^ account very well for numerous experimental situations 0. Particularly inter- 
esting is the case of a coherent state |z), a\z) = z |z), where z is a complex number. It 
corresponds to the radiation field of an ideal laser, and with Eq.(^) one easily gets Poisson 
counting distribution. 
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Pm = ^e-'\ J = v\z 
ml 



which describes a minimally bunched light source. 

In this paper we discuss measurement of electric current in the quantum regime, having 
as a primary goal a complete description of charge fiuctuations, rather than developing mea- 
surement theory for that case. We shall derive a microscopic formula for electron counting 
distribution (see Eqs.(|8D,(p!9|)) that can be used for any system, e.g., with an interaction 
or with a time- dependent potential. As an example, we find statistics in a single channel 
ideal conductor for non-equilibrium and equilibrium noise at finite temperature, and for 
zero-point equilibrium fiuctuations. 

Motivation. Similar to the photon detection, electric measurement is performed on a sys- 
tem containing an enormous number of particles — in this case fermions — and thus one 
expects the effects of Fermi statistics to be important. Also, the duration of electric mea- 
surement is typically much longer than the time it takes the system to transmit one electron 
by microscopic tunneling, scattering or diffusion. Apart from these similarities, there is, 
however, a crucial difference from the photon counting: at the detection of a current fluc- 
tuation the number of electrons is not changed, since electric charge is conserved. This has 



to be contrasted with absorption of photons in photo-counters. Related to this, there is 
another important difference: at every detection of a photon, its energy tiLV is taken from 
the radiation field, which makes plain photodetectors insensitive to zero-point fluctuations 
of electromagnetic field. On the contrary, the measurement of current fluctuation is usually 
performed without changing energy of the system, which makes the zero-point noise an un- 
avoidable component of any electric measurement . (Let us emphasize that the difference 
has nothing to do with the type of quantum statistics, Fermi or Bose, rather it is the dif- 
ference between the two kinds of measurement, e.g., see [Q, where the counting of fermions 
with a counter of "optical" kind was discussed.) 

In the classical picture, the measurement gives the charge Q{t) = fQJ(t')dt' transmitted 
during the measurement time t. The probabilities Pm of counting m electrons can then be 
obtained by averaging S{Q(t) — me) over the state of the system. In a quantum problem 
electric current is an operator, and since currents at different moments do not commute, 
the operator of transmitted charge Q{t) = jQJ{t')dt' generally does not make any sense. 
Instead, since we are interested in higher-order statistics of current fluctuations, beyond 
{j{t)) and ([j(ti)j(t2)]+), in order to compute electron counting ditribution, we have to 
include the measuring system in the quantum Hamiltonian. Our approach is motivated 
by the example of the quantum mechanical systems with strong coupling to macroscopic 
environment, introduced by Leggett, that can be treated consistently only by adding the 
"measuring environment" to the quantum problem [^. 

For that we introduce a model quantum galvanometer, a spin 1/2 that precesses in the 
magnetic field B of the current. For a classical system, the rate of precession is proportional 
to B{t), and B{t) is proportional to the current I{t): B{t) = const I{t). Therefore, the 
precession angle of the spin directly measures transmitted charge 5Q = JQl{t')dt' . We adopt 
the same measurement procedure for the quantum circuit, i.e., we include in the electron 
Hamiltonian the vector potential due to the spin: 

A{r) = -fi a X V— . (4) 

\r\ 

Thus we obtain a Hamiltonian describing motion of electrons, the measuring spin, and their 
coupling. Now, according to what has been said, we have to solve dynamics of the spin in 
the presence of the fluctuating current, find the distribution of precession angles, and then 
interpret it as a distribution of transmitted charge. Of course, a question remains about the 
back effect of the spin on the system, as in any other problem of quantum measurement. 
However, as we find below in ([13D and (|T^), only the phase of an electron state is affected 
by the presence of the spin, not the amplitude. Moreover, the phase will change only for 
the transmitted, but not for the reflected wave. As a result, the probabilities we obtain do 
not depend on the coupling constant of the spin. This justifies the assumption that the spin 
measures charge transfer in a non-invasive way. 

It is worth remarking that our scheme resembles the "Larmor clock" approach to the 
problem of traversal time for motion through a classically forbidden region |^. In this 
problem one is interested, e.g., in the time spent by a particle tunneling through a barrier. 
The Larmor clock approach involves an auxiliary constant magnetic field B added in the 
classically forbidden region, and a spin 1/2 carried by the particle that interacts with the 
field: Timt = —o'zB. The precesson angle of the spin measures traversal time. In this 



approach the field B plays a role somewhat similar to the parameter A in the generating 
function (H). The main difference from our scheme is that, in the many-particle problem, 
instead of letting each particle carry a spin, we use a single static spin coupled to current. 
Measurement scheme. Having clarified our motivation, we proceed semiphenomenologically 
and choose a new vector potential in the spin-current interaction —-jA. We replace the 
Ampere's long-range form (^) by a model vector potential 

Mr) = ^ ^^ ^^0{f{r) - /o) (5) 

concentrated on some surface S defined by the equation /(r) = /q. Here $o = hc/e, A is 
a coupling constant, /(r) is an arbitrary function, and, as usual, 9{x) = 1 for x > 0, for 
X < 0. The surface S defines a section of the conductor on which the interaction is localized: 

/■ 1 ?^ ,o Xh ^ f 

riint = / — jA d r = -—-aJs , 
J c 2e 

where 1$ = JsJds, i.e., the spin now is coupled to the total current through the section 
S. With the choice (^ of the vector potential one can study current fiuctuations in an 
arbitrary section of the conductor. Another advantage of the phenomenological Eq.(^ is 
that it involves only one Pauli matrix, which makes the spin dynamics essentially trivial. 
The choice of the quantization axis of the spin is arbitrary since (|^) will be the only spin- 
dependent part of the Hamiltonian. Finally, by switching from the smooth function (^) to 
the singular form (j^) we enforce integer values of counted charge. To understand this, let 
us note that in the "fuzzy" case (^ the measurement can start at the moment when one of 
the electrons is located somewhere in the middle of the volume where A y^ 0, and then a 
fractional part of electron charge will be counted. On the contrary, in the "sharp" case (^, 
at any time, each electron can be either on one side of the surface S or on the other side, 
and thus the charge is detected only when electron's path crosses S. We shall see below 
in a microscopic calculation, that integer values of charge follow automatically from gauge 
invariance, since the form (^) is a gradient of a scalar. 
Thus we come to the Hamiltonian 

^ ^ e " 
Ha = H{p, r), Pi=Pi Ai , (6) 

where the spin-dependent A is taken in the form (^. An essential feature of our approach 
is that we allow the constant A of coupling between the spin and the current to vary, i.e., we 
consider the spin precession as function of the parameter A. The reason is that, unlike the 
photon counting problem, our measurement scheme directly generates the function xW^ 
and then the counting probabilities Pm are obtained by reading Eq.(0) backwards. 

At this point we are able to formulate our main result. Let us construct a new Hamilto- 
nian 

-Hx = niP,r), P,=Pi- hh V,e{f{r) - /o) , (7) 

simply by supressing a^ in Eq.(^). The Hamiltonian Tix involves only quantities of the 
electron subsystem. Below we show that by studying precession of the spin coupled to the 
current, one can get the quantity 



;^(A) = (e^W-A^e-^^^* ) . (8) 

Here the brackets (...) stand for averaging over initial state of electrons. Note that xW is 
written in terms of a purely electron problem, not involving spin variables. We shall find 
that the function xW defines the result of any measurement of the spin polarization at 
the time t when the spin-current coupling is switched off. Moreover, we shall see that the 
function (|^) has the meaning of a generating function of electron counting distribution, i.e., 
Fourier transform of xW gives counting probabilities, entirely analogous to (^. 
Spin dynamics. Our goal now will be to express evolution of the spin in terms of quantities 
corresponding to the electron system. The interaction is given by Eqs.(|]),(|D. Suppose 
that the measurement starts at the moment and lasts until time t, i.e., the spin-current 
interaction is switched on during the time interval < r < t. Let us evalute the density 
matrix Ps{t) of the spin, right after it is disconnected from the circuit. We have 

Ps{t) = tnie-^'^^'pe''^^') , 

where p is initial density matrix pe®Ps at t = 0, pe is initial density matrix of electrons, and 
tre(...) means partial trace taken over electron indices, the spin indices left free. In terms of 
the spin variables, the operator e""^"^"* is a function only of ct^, and hence it is diagonal in 
spin: (T|e~*^'^*|i) = (ile^*^""*!!) = 0. In other words, if initially the spin is in a pure state, 
up or down, it will not precess. For Ps(t) this remark yields: 



Ps{t) 



^tt(o) x(A)Pn(0) 
X(-A)pit(0) pu(0) 



(9) 



Here x{^) = tr(e*''^-^*/3ee~*''^^*), where e~*^^* is the evolution operator for the problem (^. 
Now, after the spin degrees of freedom are taken care of by (||), we are left with a purely 
electron problem, that involves only electron degrees of freedom but not the spin. By using 
cyclic property of the trace one can show that x(A) in Eq.(^ is identical to (H). 

In principle, any entry of a density matrix can be measured, hence the quantity x(-^) 
is also measurable. In order to make clear the relation of x(-^) with the distribution of 
precession angles, let us recall the transformation rule for the spin 1/2 density matrix under 
rotation by an angle 9 around the 2;— axis: 



ptt ^ Pn 



(10) 



T^e{p) 
By combining this with Eq.(|D we write Ps{t) as 

Ps{t)=Y.Pmne=mx{p) , (11) 

m 

which assigns to Pm the meaning of the probability to observe precession angle m\. Let us 
finally note that such interpretation of Pm is consistent with what one expects on classical 
grounds, because for a classical magnetic moment a interacting with the current according 
to (^ the angle 6 = X corresponds to the precession due to a current pulse carrying the 
charge of one electron. 



Single- channel conductor: general formalism. In order to see Eq.(P) working, let us consider 
an ideal single channel conductor, i.e., the Schrodinger equation 

in one dimension, where the potential U{x) represents scattering region and the vector 
potential is inserted according to (H) at the x = section. In order to describe transport, we 
shall use scattering states, left and right. Their population n^^;^{E) are equilibrium Fermi 
functions with temperature T and chemical potentials shifted by eV , fJ^L — fJ'R = ^V, thus 
representing a dc voltage. We would like to study the case of small T, eV -C Ep, when only 
the states close to the Fermi level are important. Such states have almost linear dispersion 
law, and thus all wavepackets travel with the speed vp without changing shape. In such a 
case, instead of the usual scattering states, it is convenient to assume that the dispersion is 
stricktly linear and, following Landauer and Martin 0], to use the representation of time- 
dependent scattering wave packets 



i'L,T{x,t) 



6{x_), t <T 

e^^/ML(5(x_) + Bl6{x+), t > T 



where x± = x ± VF{t — r). Here r is the moment of scattering of a packet, ^l(r) and 
Bl{r) are scattering amplitudes at A = 0. To make expressions less heavy, we drop electron 
spin. Also, we assume that the scattering amplitudes are energy-independent, which is 
equivalent to replacing the barrier U{x) by Uo6{x) and is consistent with the closeness 
to Ep. The phase factors e^*'*'/^ in (|13|) are found immediately, by observing that the 
vector potential in the Schrodinger equation can be eliminated by the gauge transformation 
^(x) — > exp(zA/2 9{x))ip{x). Scattering amplitudes form a unitary matrix: 



S, 



Bl c-'^I^Ar 



(14) 



Second-quantized, electron states ( ]T3| ) lead to ip{x,t) = ipL^x^t) +ipR{x,t) with 

i^L(R){x,t) = J2'^L{R),r{x,t)c^2),T , (15) 

T 

where Ci^r and C2,r are canonical Fermi operators corresponding to the states (|T3|), the left 
and the right respectively. This defines ip{x) in (0). 

In the basis of the wavepackets (0),(0) it is straightforward to write the many-particle 
evolution operator through the single-particle scattering matrix Sx: 



^-in,t 



rt 

exp/ dTY^\ia[Sx]ijclrCj,r , (16) 



To verify (0) let us note that in the wavepacket representation (^) Fermi correlations occur 
only for the pairs of left and right states that scatter at the same instant of time. For each 
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of such pairs the evolution operator e~*^^* is 1 if both states are occupied or both are empty, 
otherwise it is given by the single-particle scattering matrix ([l^)- 
Using similar arguments, we compute 



rt 
JH-xt^-iHxt 



exp / dr J2 Wijcf^^Cj^r , (17) 



where e^ = S_lS\ is readily obtained from (|l^ : 

'c'^IAlI^ + IBlI^ 2ismXALBn 
. 2i sin XBrAl e-'^\AR\^ + \BR\^ 



e"^ 



(18) 



Using unitarity of e^ and commutation rules for Ca^r one can rewrite (|T7|) in terms of normal 
ordering: 



^i'H-xt^-i'Hxt 



exp|rfr^[e^-l],,<,c,,.: . (19) 

«J 



This form is ready to be plugged into Eq. (H) and averaged over initial state. Let us note the 
striking similarity of the two formulas obtained by different means: the fermionic Eq.(|T^ 
and the bosonic Eq.(|D. Also, let us mention that the periodicity of the matrix (|18D in A 
ensures periodicity of x(A), and thus guarantees integer values of charge. 
Single- channel conductor: results. Let us start with a simple example of a single particle 
in the state c^'T-lvac) that corresponds to scattering at the moment r. In this case, from 
(|1|) and d) one gets x(A) = e'^\A\'^ + \B\'^ for < r < t, 1 otherwise. {\A\ = \Al\ = \Ar\, 
\B\ = \Bl\ = \Bfj\) Evidently, according to Eq.(0), this simply means that for the scattering 
occurring during operation of the detector, the counting probabilities are identical to the 
one-particle scattering probabilities, as it should be expected. 

Now, let us consider current fluctuations in an equilibrium Fermi gas. Assume perfect 
transmission: -Bl(/?) = 0. Then Eq.fpTSf) gives W = iXcr^, and thus Eq.(|l^ becomes 



Jo ' ' ' ' 



i.e., the right and the left states separate. We observe that the averaging of ( poD over the 
Fermi ground state is identical to that performed in the orthogonality catastrophe calculation 

§], so we get 

X(A) = e-^^/(*'^) , (21) 

/(t,r) = (((/o*c+,ci,.rfr)2)) = I "fljl^^^ l^ ^'^. According to (D, this leads to gaus- 
sian counting statistics. 



Let us remark that, incidentally, Eq.(21) also gives a solution to another problem: the 
statistics of the number of fermions inside a segment of fixed length in one dimension. The 
relation is immediately obvious after one assigns to r in Eq. (|20D , the meaning of a coordinate 
on a line. Thus, the statistics are gaussian. 



Now, it turns out that the general case B ^ can be reduced to (pOD by a canonical 
transformation of Ca^r that makes the quadratic form in ([T7| ) diagonal. The transformation 
is related in the usual way with the eigenvectors of the matrix W. Thus, we come to 
Eqs.([20|),(pl|) with A replaced by A^,: sin ^ = 1^41 sin|. The counting statistics in this case 
are non-gaussian: 

X(A) = e-^*^(*'^) . (22) 

One checks that the second moment of the distribution agrees with the Johnson-Nyquist 
formula for the equilibrium noise. 

Finally, let us find statistics for the non-equilibrium noise. In this case, due to the 
asymmetry in the population, n^R^^E) = {exp{E±^eV)/T + 1)~^, generally one cannot 
uncouple the two channels by a canonical transformation. We calculate the statistics within 
an approximation that neglects by the effect of switching at r = and t = t. Let us 
close the axis r into a circle of length t, which amounts to restricting on periodic states: 
^(r) = ^(r ± t). For the t— periodic problem, by going to the Fourier space, one has 

X(A) = ntl + \A\\e-'^ - l)nL{Ek){l - nR^Ek)) + \A\\e'^ - l)n^(Efc)(l - n^l^fc))] , 

k 

where Ek = 27ihk/t, k is an integer. For large t, t ^ h/T or t ^ h/eV, the product is 
converted to an integral: 

ln(x(A)) = TTT / ^^ 1^(1 + l^l'(e~*' - 1)^^(1 " ^ir) + \A\\e'^ - 1)^^(1 - ul)) . 

We evaluate it analytically, and get 

x(A) = exp —tTu^u^/h , (23) 

where u± = f ±cosh~ (|y4p cosh(f + iA) + |-Bpcoshf), v = eV/2T. The answer simplifies in 
the two limits: T 3> eV and eV ^ T. In the first case we return to the equilibrium result 
(^). In the second case, corresponding to the recently discussed quantum shot noise 0, we 
have 

X(A) = {e''^\A\^ + \B\Y''^'/\ e = sgn{V) , (24) 

Analyzed according to Eq.(0), this xW leads to the binomial distribution 

rNyrn) — P Q '^n 5 

p = |Ap, q = l-Bp, N = e\V\t/h. One checks that the moments (k) = pN and ((fc^)) = pqN 
correspond directly to the Landauer formula and to the formula for the intensity of the 
quantum shot noise [^. The correction to the statistics due to the switching effects is 
insignificant [p!0|] . 

In conclusion, we introduced a quantum-mechanical scheme that gives complete statis- 
tical description of electron transport. It involves a spin 1/2 coupled to the current so that 
the spin precession measures transmitted charge. The off-diagonal part of the spin density 
matrix, as a function of the coupling constant, gives generating function for the electron 
counting statistics. We find the statistics in a single-channel ideal conductor for arbitrary 
relation between temperature and voltage. For a perfectly transmitting channel the count- 
ing distribution is gaussian, both for zero-point fluctuations and at finite temperature. At 
constant voltage and low temperature the distribution is binomial. 
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